Collective vs local measurements in qubit mixed state estimation 
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We discuss the problem of estimating a qubit mixed state. We give the optimal estimation that 
can be inferred from any given set of measurements. For collective measurements and for a large 
number, A'^, of copies, we show that the error in the estimation goes as 1/A'^. For local measurements, 
we focus on the simpler case of states lying on the equatorial plane of the Bloch sphere. We show that 
the error using plain tomography goes as l/A^^^*, while our approach leads to an error proportional 
to l/N^/\ 



I. INTRODUCTION 

Knowing the state of a system is of paramount impor- 
tance in quantum information. Quantum measurements 
provide only a partial knowledge of a state. Such a state, 
can only be reasonably reconstructed if a large number, 
A'', of identically prepared copies of the system is avail- 
able. Since the seminal work of Holevo there has been 
a lot of research on this subject. Most of the quantitative 
analysis have mainly focused onpure states, for which 
the optimal strategies IE ^ IM IS have been identi- 
fied. They give the ultimate limits that can be achieved 
in state reconstruction. However, they involve collective 
measurements that, although very interesting from the 
theoretical point of view, are very resource consuming 
and very difficult to implement in a lab. 

In the real world pure states are very scarce and so 
mixed state estimation is not just an academic issue. For 
instance, it is important to estimate the purity of a state, 
since this parameter often determines its utility to per- 
form quantum information tasks. In quantum tomog- 
raphy, a quorum of local observables is measured on a 
(large) number of copies of a state p. From the relative 
frequencies of the outcomes, one then obtains an approx- 
imation or guess pg to the signal state p. However, the 
statistical deviations often yield unphysical states, e.g. 
trpg > 1. In this case, one can either discard the results 
or use a maximum likelihood data analysis 1^. Using 
this analysis one infers the physical state that provides 
the closest theoretical probabilities to the observed fre- 
quencies. Many variants of these techniques can be found 
in the literature 0, , but there is a notorious lack of 
quantitative results (see though 

The large numbers law ensures that with an infinite 
number of copies {N — > oo) and infinite measurements 
the state could be exactly reconstructed by any sensi- 
ble method. In practice, however, one has access only 
to limited resources and the crucial issue is to quantify 
the quality of the reconstruction procedure. This is the 
question we address here. We focus on qubit states and 
use the fidelity as figure of merit. We obtain the best 
estimate for any given measurement and compute the 
analytical expressions of the average fidelity for both col- 
lective and local (von Neumann) measurements in the 



asymptotic limit (large N). 



II. BURES METRIC 

The estimation procedure goes as follows. After we 
have measured on the A^ copies of the system, some re- 
sult is obtained, which we symbolically denote by x. Note 
that X stands for both a single outcome of a collective 
measurement or a list of A^ outcomes (one for each in- 
dividual measurement) in local schemes. Based on x, 
an estimate for p can be guessed, pg{x). The fidelity is 
defined as 



/ tr \J pl'^ix) p py^{x). 



(1) 



It determines the maximum distinguishability between p 
and Pg{x) that can be achieved by any measurement |l2| . 
For qubits, Eq. 1 reads 



1 + f • R{x) + Vl - r Vl - 



(2) 



Here r = |r| and R ~ \R\, where f and R{x) are the 
Bloch vectors of the states p and pg (x) respectively [p — 
{1 + f- (j)/2; a = {cx, cTy, (7z) are the standard Pauli 
matrices] . 

The fidelity can be viewed as a "distance" between 
two density matrices. The corresponding metric is usu- 
ally known as Bures metric. From the infinitesimal "dis- 
tance" /(r, f+dr) it is easy to obtain the volume element 
(normalized to unity: J dp = 1) 



dp — — 



=dn, dn 



sin 9d9d(j) 
4^ ' 



(3) 



where dn is the invariant measure in the two-sphere. 
Eq. 13 is the natural uniform probability distribution 
function, or the a priori probability distribution for a 
completely unknown qubit state p. For our discussion 
below, we will also need dp when the density matrices 
are known to lie in a great circle of the Bloch sphere. It 
reads 



dp 



rdr 



-.de. 



(4) 
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Although in Section Hvl we use (|3J) and 10}, our main 
results in the following section are independent of any 
particular choice of the a priori distribution. 



III. FIDELITY AND OPTIMAL GUESS 

The average fidelity, hereafter fidelity in short, is the 
mean value of ^ over the a priori distribution and over 
all possible outcomes x, 



F = 



J dpf{f,R{x))p{x\f), 



(5) 



may look counterintuitive at first sight. Although the 
space of mixed states seems to be larger, they are less 
distinguishable than pure states. The fidelity Q has a 
minimum value (1 — r)/2 which is never zero but for pure 
states (r = 1). Thus the average fidelity could, in princi- 
ple, be larger than that of pure states alone. Note that 
any estimated mixed state pg has some overlap with the 
signal state p. To be more concrete, imagine one does 
random guessing, without performing any measurement 
at all, i.e., p{x\f) is uniform. Then, using Bures volume 
element the average fidelity is Fiand = l/2-|-8/(97r^), 
which is larger than the random value (F — 1/2) for pure 
states. 



where p{x\r) is the conditional probability of obtaining 
outcome x if the signal state has Bloch vector r. These 
probabilities are determined by the expectation values 
of positive operators 0{x), such that J^x^i^) — i-^-' 
p{x\r) — tr [0{x)p]. Our aim is to maximize |(SJ). 

For a given measurement scheme, {0{x)}, there al- 
ways exists an optimal guess, as we now show. We first 
introduce the four dimensional Euclidean vector 



(V 1 — r^, f). 



(6) 

+ r ■ f' and Irl = 



Note that r • r' = \/l - r^VT^ 
-y/r • r = 1. With this, the average fidelity reads 

1 + r • R(a;) 



F = Y^jdp. 



-p{x\r). 



(7) 



A straightforward use of the Schwarz inequality gives an 
upper bound of F that is saturated with the choice 



Yi{x) 



\V{x)\ 



W{x) = / dpvp{x\r). (8) 



Using the maximum fidelity is 



F = - 
2 



(9) 



Since the guess ((HJ satisfies |R(a;)| = 1 and its first com- 
ponent is non-negative, it always gives a physical state. 
In fact, for any set of measurements and any a priori dis- 
tribution, (jSJ is the best state that can be inferred and 
is the maximum fidelity. 

As the number of copies of the system becomes asymp- 
totically large, any reasonable estimation scheme leads 
to a perfect reconstruction of the state, i.e, F 1. For 
a large, but finite N, the relevant issue is knowing the 
rate at which the perfect estimation limit is attained. 
For pure states, it is well known that the best collec- 
tive strategy yields F~l-l/iV(F = l- l/(4iV) for 
states on the equator of the Bloch sphere) . It has also 
been shown recently that this asymptotic limit can be 
achieved with local measurements [g- For mixed states 
much less is known. Most of the pure state results can 
not be extrapolated to the mixed case, and some others 



IV. RESULTS 



Collective measurements 



As a first application of the results of the previous 
section, let us obtain the asymptotic behavior of the fi- 
delity with the optimal collective measurement scheme. 
The main results are contained in (9^, where an optimal 
(and minimal) generalised measurement was obtained for 
qubit density matrices and generic isotropic probability 
distributions. However no definite form for this distribu- 
tion was assumed and no explicit results were obtained. 
Our approach enables us to cast the expressions in 
in a more transparent way as well as to simplify some of 
the derivations there. The optimal measurement is repre- 
sented by a set of positive operators and conditional prob- 
abilities that can be labelled by two indices x — {k,rn). 
The discrete index k labels the representations of the 
symmetric space spanned by {p®^} onto which the pos- 
itive operators of the measurement project, whereas the 
unit vector rn labels a continuous set of outcomes in the 
two-sphere [I^- We have 



p{k, m\r} = Cfe 



1 



N/2-k 



1 + f* • m 



2k 



where 



Ck 



f N \ {2k + 1)^ 
yN/2 + k ) N/2 + k + l' 



From © and ^ we obtain 



F 



dm |V(fc, m) 



(10) 



(11) 



(12) 



The sum in (|12(l runs from = (fc = 1/2) for even 
(odd) to fc = N/2. Taking advantage of the rotational 
invariance, the integrals in (|12|l can be evaluated ex- 
actly. The computation of the asymptotic limit is rather 
lengthy and will not be reproduced here The final 
result is 



F = 1 - 



3 

'SttJ N 



(13) 
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Note that this fidehty is only shghtly worse than that of 
pure states: 3/4 + 4/(37r) = 1.17 > 1. 

The expression H13(l also gives us important informa- 
tion about the optimal fidelity when the a priori proba- 
bility distribution corresponds to states known to lie in 
the equator plane of the Bloch sphere Q). Since in this 
situation we have more information about the states, the 
fidelity cannot be worse than H13|l . i.e. the error, defined 
a,s E — 1 — F, must satisfy E < ^/N, where f is a con- 
stant. 



B. Local measurements 

Let us now tackle the problem of reconstructing a qubit 
state p from local measurements. The analytical ex- 
pressions turn out te be rather difficult to obtain and 
for simplicity only the case of states that are known to 
lie on the equator plane of the Bloch sphere ^ will be 
considered in this note. This is a non trivial case that 
can be relevant for quantum optics (e.g., for polarization 
states of photons). We will only sketch our main results. 
The techniques we have used are essentially contained 
in 1^ . Full details of the calculations will be presented 
elsewhere 14] . 

Consider N — 2Jsf copies of the state p. Quantum 
state tomography tells us that von Neumann measure- 
ments along two fixed orthogonal directions, x and y, 
are sufficient to reconstruct the state. After the mea- 
surements, we obtain a set of outcomes -1-1 and —1 with 
relative frequencies a; and 1 — cti, respectively {i — x.y). 
This occurs with probability 



p{a\r) 



n 



1 



7^ a. 



1 



■ (14) 



In quantum tomography the guess is given by 



Rt(<3) = (Vl — i?^, i? COS7, i? sin7), 
i?cos7 — 2ax — 1, i?sin7 = 2ay — 1. 



(15) 



In many instances, however, the statistical fluctuations 
produce an unphysical guess (the square root term be- 
comes imaginary). If one discards these cases, the asymp- 
totic behavior of the fidelity can be shown to be: F — 
1 - Ct/^^/" + where is a constant. Although 
plain tomography yields a perfect reconstruction of the 
state in the asymptotic limit, it is much worse than the 
optimal collective scheme (note the power 1/4 of iV as 
compared to the power 1 in Eq. I13|) . One may suspect 
that the culprit of this behavior is the number of copies 
discarded, but we now show that it is not entirely so. 

Within the maximum likelihood framework 7] all 
available data is used. If i? < 1, the guess is the to- 
mographic one: Rml = R-t (see Eg. I15|l . and 



if i? > 1, where $ is the solution of the equation 
cos(2$) = i?cos(7 -I- <&). In the asymptotic limit one 
can expand this equation as a power series in [R — 1), 
<I> = 7— (i?— 1) cot 27 -(-••• . In fact only the first term is 
necessary for our calculation. After some effort one gets 



F = 1 - 



Ar3/4 



ML 



£(1/4)3 

25/497r2 



0.2256. (17) 



Notice the significant increase in the rate at which the 
fidelity approaches unity as compared to plain tomogra- 
phy. 

Finally, we have computed the fidelity for the optimal 
guess 0. Here again, all available data is used to pro- 
duce a reconstruction of the state. In Fig. 1 we compare 
the optimal guess and maximum likelihood methods for 
up to = 20 copies of a state. It is clear that the opti- 
mal guess strategy always performs better. We have also 
obtained the asymptotic limit. The fidelity reads in this 



F = 1- 



7V3/4 



Co - 0.1708, 



(18) 



Rml = (O,cos($),sin($)) 



(16) 



where can be computed analytically (see appendix). 
Notice that this fidelity approaches unity at a rate sim- 
ilar to the maximum likelihood one (|17|l . but the coef- 
ficient of the first correction is lower (^o < Cml), as 
it should. The most important parameter is the expo- 
nent of the 1/A^ term in H18() . It shows that, there is a 
gap in the quality of the reconstruction process between 
fixed local measurements and optimal collective schemes 
(recall Eq. I13|) . One may argue that we have not ex- 
ploited classical communication, i.e, we have not designed 
each individual measurement according to the outcomes 
of the previous ones. We have also explored this pos- 
sibility numerically and observed that the fidelity is al- 
most identical to that obtained from the optimal guess 
with measurements along two fixed orthogonal directions. 
Therefore, we are led to conjecture that an error rate 
E r^l /A^3/4 ^.j^g lowest that can be achieved using any 
local scheme. 



V. CONCLUSIONS 

We have obtained the optimal reconstruction of a gen- 
eral qubit state for any given set of measurements and 
have illustrated our results in some interesting cases. We 
have computed the asymptotic expression of the fidelity 
for the optimal collective scheme. For local measure- 
ments we have considered the simpler but important case 
of states lying on the equator plane of the Bloch sphere. 
We have shown that the performance of plain tomogra- 
phy is very poor, with an error that goes as E ^ X/N^/'^ 
for large N . We have shown that maximum likelihood 
does provide a much better estimation: E ^ Us- 
ing the same data, the optimal guess analysis gives the 
best reconstruction of the signal state. Despite this im- 
provement, the asymptotic behavior of the fidelity does 
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FIG. 1: Average fidelities in terms of the number of copies. 
The diamonds (crosses) correspond to the maximum likeli- 
hood (optimal) guess. 
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APPENDIX 



The constant in H18|l is the sum of three terms that 
can be written as 



not saturate the optimal collective bound. This is in con- 
trast to pure state estimation, where local measurements 
can perform optimally in the asymptotic regime. Al- 
though we have mainly focussed on measurements along 
fixed orthogonal directions, we have also analysed the 
most general local strategy, in which one is entitled to 
change these directions after each individual measure- 
ment. Our results strongly suggest that for mixed states 
the asymptotic behavior of the optimal collective schemes 
cannot be attained by any local strategy. 
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io = -\i^{Ah+b,-V2b3); 



hi 



487r 

OO 

dx 



V2 



\x^/HuY[Ki/i{-x)] ^x^/^ 
e^[erfc(^/2^) - 4]erfc(V2x) 



dx 



dx 



a;3/4lm[ifi/4(-a;)] 
e^erfc^(\/2lF) 



/o xViK,„{x) ' 
where K^, and erfc are the modified Bessel function and 
complementary error function respectively |l5j| . The in- 
tegrals can be evaluated numerically: hi = 0.197241, 
&2 = 1.61451 and 63 = 0.31400. With these numbers 
we have ~ 0.17083. 
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